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Abstract 


The infrared absorption spectra of thin crystalline films of sulfur 
dioxide at 90 K are reported in the 2700 ott* - 450 cm"! region. The 
observed multiplicity of the spectral features in the regions of fundamen- 
tals is attributed to factor group splittings of the modes in a biaxial 
crystal lattice and the naturally present minor 36$ anc j 18g isotopic 
species. Complex refractive indices determined by an iterative 
Kramers-Kronig analysis of the extinction data, and absolute band strengths 
derived from them, are also reported in this region. 


Introduction 


Sulfur dioxide has been investigated spectroscopically in exhaustive 
detail. 1- 7 Because of its importance as a pollutant in factory and power 
plant smoke, and in automobile exhausts, methods have been devised for 
remote sensing of this molecule. 2 *^ These methods are based on the charac- 
teristic absorption/emission peaks in the ultraviolet and infrared 

regions. ^ 

Such remote sensing methods have been used to identify sulfur dioxide 
on Jupiter's moon, Io. The characteristic emission at 1360 cm -1 recorded 
by the Voyager IRIS experiment was responsible for the unambiguous iden- 
tification of this gas in the satellite's atmosphere by Pearl, et al. 11 
The 2456.2 cnr 1 band observed in the reflectance spectrum of Io, recorded 
by Cruikshank, et al., 12 has been assigned to sulfur dioxide frost 12 on the 
satellite's surface. This is corroborated by a relatively broad feature, 
characteristic of solid material, in the thermal emission spectrum of Io 
recorded by the Voyager infrared experiment; positioned at 552 cm -1 ; this 
band has has been identified with solid SO 2 on the surface. 14 

Laboratory data on the extinction coefficient of the 1360 cm -1 band 
of gaseous sulfur dioxide were utilized by Pearl, et al. 11 to estimate the 
gas abundance in Io's atmosphere. Corresponding modeling of the surface 
requires laboratory data on the complex refractive indices of appropriate 
solid materials in the thermal infrared region. No such data for solid 
sulfur dioxide are available in the literature, although extensive investi- 
gations of the vibrational infrared and Raman spectra of crystalline 
sulfur dioxide have been published. i* 2 *^, 7, 15, 16 

In this report, we present a detailed interpretation of the infrared 
active modes of crystalline sulfur dioxide, thereby clarifying a number of 


inconsistencies in previously reported work. From quantitative 
transmission data, the complex refractive indices and the integrated 
extinction coefficients of the absorption bands are also obtained. 

Experimental 

The apparatus and procedure are similar to those described in an 
earlier report on crystalline 0 , 2 ^ 2 .'^ Sulfur dioxide obtained from 
Matheson (99.5%) was further purified by vacuum distillation to remove the 
main impurity, carbon dioxide. Thin films (0.27-2.8 y thickness) were 
deposited on a KRS-5 substrate at 90 K and annealed for 30 minutes. The 
thicknesses of the films were calculated from the channel fringes in the 
near infrared, using a refractive index of 1.53 (based on that of the 
liquid, adjusted by using the Lorentz-Lorenz formula). 1 ** Quantitative 
infrared transmission spectra were recorded on a Perkin Elmer 1800 FTIR 
instrument. A spectral resolution of 0.6 cm " 1 was employed throughout the 
measured 2700-450 cm -1 spectral range. Characteristic spectra of the thin 
film samples are reproduced in figure 1 . 

Data Analysis 

An isolated SO 2 molecule (point group C 2 V ) has three fundamental 
vibrational modes. In the gas phase the frequencies and the classification 


(A 2 ) : 1151 cm ' 1 
$2 (Aj) : 519 cm -1 
$3 (Bj) : 1361 cm * 1 


of these modes are:** 


In the solid phase the intermolecular interactions produce significant 
changes in positions and strength of the bands relative to those in the gas 
phase. Bands due to the naturally present isotopic species -^S02 (4.4%), 
34$18oo (0.41%) and ^®S02 (0.07%) are also observed in both phases. 

The crystal structure of SO 2 is orthorhombic (space group C^y, Z=2) 
with the molecules occupying C 2 sites. ^ The correlation diagram (figure 
2) shows the expected multiplcity of the molecular modes due to crystal 
field effects. Thus, for 01 and 9 2 modes only one of the two crystal 
field split components (Aj) is infrared active whereas 03 is split into 
two active components (Bj and B 2 ). Because of the biaxial nature of the 
crystal each infrared active mode is further split into 2 components, TOj 
and TO 2 . 2 O This is precisely what is observed in the spectra. 

Peak frequencies derived from our measurements are listed and com- 
pared with earlier results in Table 1. In the regions of fundamentals 
our observed frequencies are very similar to those reported previously. 
Earlier works have attempted to explain the infrared bands in terms of 
the normal modes of SO 2 and it isotopic species. Brooker* 5 discussed the 
Raman spectrum of crystalline sulfur dioxide in terms of Transverse (TO) 
and Longitudinal (L0) modes but failed to recognize the splitting of the 
TO components in a biaxial crystal. The assignments proposed here are 
consistent with the predictions based on the crystal field effects, 
the expected frequency shifts and relative intensities of the naturally 
present isotopic species. A In the following discussion the spectrum is 
divided into four regions. 

a) 510 - 540 cm" 1 

The strongest bands, centered at 520.9 cm"* and 522.7 cm - *, 
assigned to the 2 TO components of the $2 ( A j) mode. The two weak 


are 



shoulders peaked at 516.7 cm“l and 514.0 cm"* are associated with the $2 
modes of S* 8 00 and 3 ®S02 species respectively. The observed shifts of 
these peaks from the average frequency of the two main TO components (5.1 
cnT* for S^OO and 7.8 cm"* for 38 S02) compare favorably with the 
corresponding calculated values^ of 4.6 cm - * and 10.6 cm - *. 

b) 1120 - 1150 cnT 1 

Results in this region follow from arguments similar to those for the 

$2 region because the and $2 modes have the same symmetry characters. 

The two TO components are associated with the 1140.0 cm -1 and the 1143.0 

cm - l peaks; in all previous work the band at 1140.0 cm _ l had been assigned 

to of 34$02. However, its observed intensity relative to that of the 

2 ( SO2 ) is much larger than the expected 4.4%; in addition, the shift of 3 

cm~l from the main peak at 1143.0 cnT* is substantially smaller than the 

calculated value, 7.3 cm -1 . Here we associate the weak shoulder at 1134.0 

cm -1 with 34 S02; we associate the shoulder at 1120.5 cm -1 with S 18 00. Both 

of these assignments are consistent with their relative intensities and the 

2 

respective calculated frequency shifts 7.3 cm -1 , and 28.2 cm -1 . 

c) 1280 - 1350 cm" 1 

As discussed above, 4 TO components (2Bj and 2B2) of 03 are expected 
in the infrared spectrum. These are, in fact, observed at 1302.8 cnr*, 
1304.3 cm”l, 1310.1 cm”*, and 1323.3 cm"l. Their classification into Bj 
and B2 species as given in the table 1 is somewhat arbitrary in the absence 
of detailed polarization data. Our results were obtained for thin films of 
randomly oriented crystals. Brooker carried out polarized Raman studies;^ 8 
however, because of partly oriented crystals and significant polycrystalli- 
nity of the sample, there was considerable spillover of modes of other spe- 


cies in a scattering geometry designed to obtain modes of a specific 
symmetry. For example, in the scattering geometry, a(cc)b in which only Aj 
modes should appear, the band at 1310.6 cm -1 assigned to 93(61, TO) appears 
strongly. We assign the observed infrared peaks at 1295.6 cm" 1 , 1291.2 
cm' 1 , and 1287.7 cm' 1 to 03 modes of 34 S02, S 13 00 and 36 S02, respectively. 
Their relative intensities are roughly in the ratios of their 
natural abundance and their frequency shifts from the average of the four 
main components (~17.5 cm" 1 , 21.9 cm" 1 and 25.3 cm -1 , respectively) are in 
agreement with the corresponding calculated values 16.9 cm -1 , 21.0 cm" 1 and 
25.4 cm" 1 . 4 

d) 2270 cm" 1 - 2460 cm" 1 

Sharp features at 2287.4 cm" 1 and 2273.9 cm" 1 are assigned to the 2b 1 
overtones of SO2 and 34 S02 species, respectively. The features at 2456.2 
cm" 1 and 2433.7 cm" 1 have proper relative intensities and spectral shifts 
to be assigned to the bj + b3 combinations of SO2 and 34 $02 species. These 
assignments are in agreement with those of Barbe, et al. 3 

In the regions of all the fundamentals, broad absorptions are observed 
on the high frequency side of the TO peaks. These are characteristic of 
strongly polar modes for which the L0-T0 splittings are large; the L0 modes 
are made infrared active by anharmonic lattice interactions. 20 In prin- 
ciple, appropriate scattering geometry can be employed to obtain L0 fre- 
quencies in the Raman spectra of polar crystals. The reported Raman 
spectra of S02» 1 ^ though not obtained for perfectly oriented crystals, give 
indications'of peaks around 1150 cm" 1 , 530 cm" 1 , and 1340 cm" 1 which may be 



associated with LO components of $2 and ^3» respectively. These fre- 
quencies are in approximately the same range as those of the shoulders in 
our infrared spectra. 

The frequencies of the A2 components of and $2 ap e not established; 
therefore, the magnitude of crystal field splittings of these modes is not 
known. As mentioned above, the classification of Bj and B2 crystal field 
split components of $3 is also somewhat arbitrary because of lack of 
polarization data. The basis, for our choice of assignments of S3 com- 
ponents (Table 1 ) is the large splitting (Bj-B 2 separation) expected due to 
dipole-dipole interaction between strongly polar modes. Strong dipolar 
interactions are also indicated by the fact that the observed LO-TO split- 
tings (~5 cm“l for Sj and $2 and ~25 cm"l for S3) are roughly proportional 
to their intensities.^ 

Optical Constants and Band Strengths 

To properly evaluate integrated strengths of the bands requires 
knowledge of the wavenumber dependent absorption coefficient. Thin film 
measurements do not provide this function directly because the transmitted 
signal is reduced by reflection losses at the sample surfaces, as well as 
by absorption within the sample. Both of these effects are related to the 
complex refractive index of the material m($) = n($) - ik($). 

To derive m(o) from the data, an iterative procedure similar to that 
suggested by Warren^* was adopted. Taking account of both reflection and 
absorption, the intensity transmitted through a plane absorbing film of 
thickness d deposited on a nonabsorbing substrate is given approximately 

by21,22 



Here 


(l-RxHl-RgHl-F^) e ^ d 
l-RiR 2 


a = 4 tt vk 

Rj = reflection coefficient at the vacuum-film interface 

R 2 = reflection coefficient at the film-substrate interface 

R 3 = reflection coefficient at the substrate-vacuum interface 


The transmittance, T, is defined as the ratio of the measured intensity to 
the baseline: 


T = 


IT 

^Baseline 


„ (1-Ri)(1-R 2 ) 

l-R r R 2 


!- r 2* r 3 e “ a True d 
(1-Rl*)(l-R 2 *) 


where * refers to conditions just outside the regions of absorption bands. 

To initiate the iterative procedure, an estimate of the true absorption 
coefficient, ajrue* obtained from the ratio of transmission spectra for 
two different sample thicknesses (the reflection losses are thereby can- 
celled, to a first approximation): 

“True * oTpjJ ,n ^ 


The corresponding k values are then used in a Kramers-Kronig analysis to eva- 
luate the real part, n, of the refractive index of the film. A new estimate of 
“True 1S then obtained by using the n and k values in the expression 

-InT _ + I In ( 1 ~ r 1*( 1 - r 2*)( 1 - r 1 r 2) 

d “True d (1-R 1 )(1-R 2 )(1-R 2 *R 3 ) 



Obtaining n through the Kramers -Kronig analysis and reevaluating ctj rue defines 
the procedure which was iterated to convergence. 

In practice, the method was only applied within absorption bands, for 
in addition to the absorption signatures of SOg, the actual data contain: 
a) weak absorption features associated with impurities within the sample 
and/or the instrument windows (principally H 2 O and vacuum pump oil, respec- 
tively; these features do not overlap those due to SO 2 ); b) a slight slope 
due to instrumental adjustment (this is not caused by scattering in the 
samples, all of which were clear and transparent to the eye; in addition, 
the signal decreases with decreasing wavenumber whereas that due to scat- 
tering would increase with decreasing wavenumber); and c) channel fringes 
due to coherent addition of radiation multiply reflected between the 
parallel surfaces of the films (these typically had periodicities of 10^ - 
10* cm"*). For the present analysis, we assume that T * 1 outside of the 
absorption bands, and define a baseline within the bands by a linear func- 
tion which is tangent to the channel fringes at the band edges. Plots of n 
and k are shown in Figure 3. As a check, the complex refractive indices 
obtained in this way were used to calculate transmission spectra for 
samples of known thicknesses; Figure 4 shows observed and calculated 
transmittance in the 1370 - 1270 cm"l region for three different 
thicknesses. In all cases the agreement with the laboratory data is within 
10%. Because of the choice of baseline, the complex refractive indices in 
the wings of the absorption bands are somewhat distorted. We consider the 
k values accurate to +_ 10% within the central regions of the fundamentals 
and to +_ 20%. for the weak 29j band. Table 2 gives selected n and k values 
in the regions of the fundamentals. 

The derived absorption coefficients were used to determine the 


integrated band strengths 


A ■ I “True<'» d0 


Results are given in Table 3, along with approximate values derived by a least 
squares fit to the transmittance data uncorrected for reflection losses: 

A* = In T(b)db 

Also included in the table are the integrated band strengths for the gas. 

The greatest source of error in these determinations is in the thickness 
measurements of the films. We estimate this error to be _+ 10%. For weaker 
combination and overtone bands there is an additional error, possibly up to 
10% in the measurement of their integrated absorbance areas. Overall we consider 
our results good to within 10% for the fundamentals and to within 20% for 
the weaker overtones and combination bands. The strengths of the $2 and ^3 
bands are not significantly changed on condensation. However, the changes 
in the strengths of the bj fundamental and the 2bj overtone are substantial, 
and lie outside the limits of experimental error. The cause of this beha- 
vior is unknown, as a clear understanding of the effect of crystal field on 
the band extinctions is not available at present. 
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TABLE 1 

Vibrational Frequencies for Crystalline Sulfur 


Dioxide; Assignments are based on C2 V Space Group 


Anderson and 
Savoie 
1965 

IR, 77°K 
cm - * 

Giguere and Falk 
1956 

IR, 98°K 
cm - * 

Barbe et al . 
1971 

IR, 77°K 
cm - * 

Brooker 

1984 

Raman, 77° 
cnT* 

This Work 
1986 

K IR, 90°K 

cm - * Assignment 



514.3 


514.0 

02(Ai)[Sl8 0 0] 


517 


517.2 

516.7 

MAi)[ 34 S02] 


521 

522 

521.9 

520.9 

522.7 

Ai(TOi) 

V2(S02) 

Ai(T0 2 ) 

524 

528 

535 


524.3 

~ 530 
(broad 

M L °) 

shoulder) 

542 






1122 

1121 

1120.6 

1120.5 

1120.5 

^>l(A 1 )[S 18 00D 



1134.9 

1134.6 

1134.0 

^i(Ai)C 34 so 2 ] 

1141 

1144 

1140 

1144 

1140.4 

1144.8 

1140.5 

1143.8 

1140.0 

1143.0 

Ai(TOi) 

$]S0 2 

Ai(T0 2 ) 

1148 

1153 



1148.2 

-1150 

(broad 

Oi(L0) 

shoulder) 



1288 


1287.7 

9 3 [ 36 S0 2 ] 





1291.2 

^ 3 [Sl800] 





1295.6 

$ 3 [ 34 so 2 ] 

1304 

1303 

1302.8 

1304 

1302.8 

1304.3 

Bi(TOi) 

^3 S0 2 

Bi(T0 2 ) 


1305.0 



TABLE 1 (Continued) 

Vibrational Frequencies for Crystalline Sulfur 


Dioxide; Assignments are based on C2 V Space Group 


Anderson and 
Savoie 
1965 

IR, 77°K 
cm" 1 

Giguere and Falk 
1956 

IR, 98°K 
cm* 1 

Barbe et al. 
1971 

IR, 77°K 
cm"! 

Brooker 

1984 

Raman, 77°K 
cm* 1 

This Work 
1986 

IR, 98°K 
cm* 1 

As si gnment 

1312 

1310 

1315.2 

1310.6 

1310.1 

B 2 (T0i) 






$ 3 S0 2 

1324 

1322 

1323.2 

1323.6 

1323.3 

B2(T0 2 ) 

1341 

1334 


1337 

'1345 

9 3 (L0) 





(broad shoulder) 

1351 



1353 





2274.3 


2273.9 

2bi( 34 S0 2 ) 


2287 

2288 


2287.4 

2^i(S0 2 ) 


2432 

2433.5 


2433.7 

$i+$3( 34 S02) 


2455 

2456.5 


2456.2 

^l+b 3 (S0 2 ) 


TABLE 2 


Selected Values of n and k in the Infrared Absorption 
Regions of Crystalline SO 2 


$ (cnf-1) 

n 1 

k 

2458.2 

1.506* 

0.025 

2456.2 

1.534 

0.062* 

2454.3 

1.559* 

0.027 

1400.2 

1.361 

0.021 

1375.0 

1.299 

0.047 

1344.7 

0.999* 

0.317 

1334.5 

1.043 

0.570 

1323.4 

1.359 

1.399* 

1316.5 

1.650 

0.842 

1310.2 

2.303 

1.410* 

1305.7 

1.799* 

0.539 

1304.5 

2.142 

0.921* 

1302.1 

2.645* 

0.370 

1210.0 

1.597 

0.000 

1151.2 

1.101* 

0.311 

1150.0 

1.110 

0.424 

1145.5 

1.302 

0.902 

1143.0 

2.178 

1.737* 

1142.5 

3.040* 

0.737 

1140.0 

2.138 

0.976* 

1800.0 

1.564 

0.00 

550.5 

1.381 

0.023 

537.3 

1.161* 

0.237 

522.7 

1.772 

1.234* 

520.9 

2.769 

1.660* 

520.2 

3.253* 

0.574 

516.7 

2.066 

0.079* 


1 Values of n are based on an assumed value of 1.53 at 2700 cm'* in the 
Kramers-Kronig analysis of the absorption data. 

Data for n and k at 0.3 citT* interval in the absorption regions are 
available from the authors. 


* Peak positions 




TABLE 3 


Integrated Band Strengths for Gaseous and Solid 
SO 2 (All Values in Units of lCr cm~2) 


Band 

A Gas 

A Sol id 

★ 

A Sol id 

*1 

7.6 

13.6 

15.2 

*2 

9.4 

8.7 

12.0 

$3 

67.8 

61.8 

68.5 

2*1 

0.79 

0.16 

0.13 

*1+S3 

0.97 

1.09 

0.97 



Figure 1. Infrared transmission spectra of thin films of crystalline SO 2 at 
90 K 

(a) 2.8 micron thickness 

(b) 0.27 micron thickness 

Transmission scales have been arbitrarily shifted for the two 
spectra. 


Figure 2. Correlation of vibrational modes of SO 2 molecules in crystalline SO 2 . 
Figure 3. Complex refractive indices of crystalline SO 2 at 90 K. 


Figure 4. Comparison of calculated (•••) and laboratory (- 
spectra in the $3 region. 

(a) 0.27 y sample thickness 

(b) 0.52 y sample thickness 

(c) 0.75 y sample 


-) transmission 


thickness 



TRANSMISSION 



WAVENUMBER 



Molecular 

(C 2 v) 


Group 


Site Group 

(C 2 ) 


Factor Group 

(C 2v ) 


Selection Rules 


i ,t>2 A 1 


A 



Al (T Z ) 

A 2 


IR , R 
(T0 1 ,T0 2 ) 
> R 


*3 B 1 



Bl T( X ) 
B 2 (Ty) 


IR , R 
(T0i,T0 2 ) 
IR , R 
(T0i,T0 2 ) 


IR: Infrared Active 

R: Raman Active 

T X , Ty, T i indicate directions of dipole change during vibration 

T0i,T0 2 : Two possible phonon propagation directions perpendicular to the 

direction of dipole change 



z 





TRANSMITTANCE 




